' AD-A033  915 


UNCLASSIFIED 


NAVAL  RESEARCH  LAB  WASHINGTON  0 C p/g  20/9 

THE  INITIAL  BEHAVIOR  OF  PLASMONS  IN  THE  PRESENCE  OF  A LONG  WAVE— ETC (U) 
NOV  76  W E HOBBS » W M MANhEIMER 

NRL-MR-3406  Ml 


l(rl 

A0AO339I5 

m 

| 

END 

DATE 

FILMED 

2 -77 

ADA033915 


rhe  Initial  Behavior  of  Plasmons  in  the  Presence  of 
a Long  Wavelength  Large  Amplitude  Phonon 


W.  E.  Hobbs,  Jr.  and  W.  M.  Manheimer 

Plasma  Dynamics  Branch 
Plasma  Physics  Division 


November  1976  ^ 


NAVAL  RESEARCH  LABORATORY 
Washington,  D.C. 


Approved  lor  public  release:  distribution  unlimited. 


SECURITY  CLASSIFICATION  OF  ThiS  RAGE  fWhon  Dote  Enlorod) 

. REPORT  DOCUMENTATION  PAGE 


la;  GOVT  ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
1.  RECIPIENT'S  CATALOG  NUMBER 


«t  Title  CawV  SyfltUfl  - . 

*EBE  INITIAL  BEHAVIOR  OF  PLASMONS  IN  THE 
fRESENCE  OF  A LONG  WAVELENGTH  LARGE 
AMPLITUDE  PHONON  ' 


-A.  I TEE  OF  R^gORT  A PERIOO  COVENEO 

Interim  repS  on  a continuing 
NRL  problem.*  * 

*.  PERFORMING  ONO.  REPORT  NUMSEP 


7 AU  ThOH'IJ 


fw.E.  Hobbs,  Jr.,*  atari  W M.Manheimer 


».  CONTRACT  OR  GRANT  NUMBER!’.,! 


» PERFORMING  ORGANIZATION  NAME  ANO  AOORESt 

Naval  Research  Laboratory  i^ 

Washington,  D.C.  20375 

It  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Energy  Research  and  Development  Administration 
Washington,  D.C.  20545 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  • WORK  UNIT  NUMBERS 


NRL  Froblerps^HO 
AEC  AT-i&IG 

U-  REPORT  OATB  ' 

Novenbar  1976 


'02-37  & H02-29B 


rj-  NUMBER  OP-WR® 

is  7r.  7 


MONITORING  AGENCY  NAME  A AOORESS (It  dllloront  from  Controlling  Oil tco)  [ 15.  SECURITY 


UNCLASSIFIED 


ISA.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


II*.  distribution  STATEMENT  tal  mi.  Report) 


Approved  for  public  release;  distribution  unlimited. 


’T  DISTRIBUTION  STATEMENT  fo I m.  MiIiki  . 


1 In  Block  TO.  II  dlllpri II  tram  twill  I K.'.-JJ.IIH 


ms  wit.:,  ’its**  sy 

ti"3  tail  SscitM  AqS 

WilMOWta  O 


i*  supplementary  notes 

*NRC  Associate.  

IT  

niiviiisrioa/miuaiiin  cma 

1*.  KEY  WORDS  rContlnu*  on  nmH  .Id.  II  AM.i.Arr  «nd  IdwilKr  Br  Slack  mmbtr)  1 C <1.  IHIt.  m^w  ^ 

Nonlinear  waves  | A 

Fluid  simulation  j / L 

Langmuir  solitons  * [f 

Solitons  4 L — . n 

20.  ABSTRACT  ( Continue  on  reeeree  aide  II  neceeomry  md  Idonlltr  by  block  number) 

•—  ^>(Jsing  a one  dimensional  two-fluid  model  the  behavior  of  electrostatic  electron  plasma  oscilla- 
tions in  the  presence  of  a large  amplitude  long  wavelength  ion  acoustic  wave  has  been  studied.  The 
electron  waves  are  observed  to  collect  initially  in  regions  of  maximum  density  gradients.  The  initial 
increase  in  the  variance  of  the  plasmon  wavenumbers  is  quadratic  in  time.  When  the  high  frequency 
fields  are  characterized  by  a narrow  spectrum  in  k-space,  a solution  of  the  fluid  equations  is  obtained 
which  exhibits  these  results. 


1 - — ■ - . 

DD  i 1473  COITION  OF  I NOV  *s  IS  OBSOLETE  1 *’?  V 7 (Aj  <"V') 

S/N  0102-014-  4601  - --fa  J ' 

SECURITY  CLASSIFICATION  OF  THIS  PAOE  fWNMi  0.1.  RnlPnd) 

f 


■' 


THE  INITIAL  BEHAVIOR  OF  PLASMONS  IN  THE  PRESENCE  OF 
A LONG  WAVELENGTH  LARGE  AMPLITUDE  PHONON 


During  the  last  few  years  there  has  been  extensive  examination 
of  the  nonlinear  system  of  Langmuir  oscillations.  Motivation  for 
these  studies  is  their  relevance  to  plasma  heating  via  powerful 
relativistic  electron  beams  or  powerful  coherent  laser  radiation. 

A particularly  striking  effect  is  the  observation  in  computer  simu- 
lations of  the  localization  of  electrostatic  fields  accompanied  by 
self-consistent  density  cavities.  Generally,  the  theories  of  this 
phenomenon  have  involved  the  nonlinear  Schrodlnger  equation  and  its 
envelope  soliton  solutions. 

To  complement  these  studies  we  have  made  numerical  studies  of 
the  two-fluid  plasma  equations.  We  initialize  the  electron  fluid 
with  a set  of  long  wavelength  electrostatic  waves  (kg  < k^/3,  where 
kfi  is  the  Initial  plasmon  wavenumber  and  k^  is  the  Debye  wavenumber). 
Inhomogeneities  in  this  field  will  lead  to  the  development  of  envelope 
Langmuir  solltons.  The  following  heuristic  description  of  the  process 

is  from  a note  by  Nishikawa  et  al.1 

Suppose  a small  amplitude  modulation  is  Introduced  to 
the  uniform  plasma  oscillations.  Then  some  of  the 
plasma  density  is  expelled  from  the  larger  amplitude 
region  by  the  pondermotive  force.  As  a result,  the 
local  plasma  density,  and  hence  the  local  plasma 
frequency  is  decreased  in  the  larger  amplitude  region. 

A decrease  of  plasma  frequency  implies  an  increase  of 
the  "wavenumber",  which  means  a further  localization 
of  the  wave  field.  In  this  way,  the  small-amplitude 
modulation  is  amplified. 

The  objective  of  our  calculations  is  to  numerically  study  this 

process.  To  get  our  system  going,  we  Initialize  the  ion  fluid  with 
Note:  Manuscript  submitted  October  29,  1976. 
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a long  wavelength  sinusoidal  Ion  acoustic  wave  (k  **  k / 20,  where  k 
is  the  wavenumber  of  the  acoustic  wave).  Our  computations  show  a 
localization  of  the  Langmuir  radiation.  However , In  contrast  to  the 
above  description,  the  observed  decrease  In  the  wavenumber  is  Initially 
localized  In  regions  of  strong  plasma  density  gradients.  We  also 
measure  the  variance  in  the  Fourier  wavenumbers  during  the  localization 
process.  This  variance  Initially  proceeds  quadratlcally  In  time,  not 
linearly  which  would  be  the  case  If  the  wave  amplitudes  were  following 
a diffusion  equation  in  k-space. 

These  results  are  understood  by  noting  that  in  our  computations 
we  are  solving  an  Initial  value  problem  and  that  for  our  Initial 
condition  the  high  frequency  fields  have  a negligible  group  velocity. 
For  this  condition  we  find  that  the  fluid  equations  may  be  solved 
analytically,  and  that  the  solution  exhibits  the  above  characteristics. 

As  the  Fourier  wavenumbers  spread,  the  amount  of  frequency 
variation  In  the  spectrum  of  the  high  frequency  waves  Increases.  In 
our  simulations  we  then  observe  the  localization  to  shift  to  regions 
of  low  plasma  density  and  the  development  of  envelope  wave  structures. 
Initially  the  effect  of  the  density  modulation  is  a dynamic  process 
of  wave  localization.  Later  the  system  evolves  so  that  the  density 
modulation  is  balanced  by  the  dispersion  or  pressure  associated  with 
a wave  amplitude  modulation. 

We  will  first  describe  our  numerical  model.  We  will  then  relate 
our  results  and  discuss  them  In  terms  of  the  fluid  equations. 

Our  model  consists  of  the  one-dlmenslcnal  fluid  equations  for 
the  electrons  and  Ions, 
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and  Poisson '8  equation, 
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The  independent  variables  x and  t have  been  scaled  to  the  inverse 
electron  Debye  wavenumber  k^1  and  the  inverse  electron  plasma  frequency 
u^1  respectively.  The  subscripts  e and  i refer  to  electron  and  ion 
properties  respectively;  m is  mass,  T is  temperature,  and  y is  the 
ratio  of  specific  heats.  The  densities,  ng  and  n^,  are  scaled  to  the 
equilibrium  density  nQ;  the  velocities,  vfi  and  v^,  are  scaled  to  the 
electron  thermal  velocity  and  the  ion  acoustic  velocity 

vv  v/nT7m^  respectively;  and  the  potential  4 is  scaled  to  the 
electron  temperature  in  volts.  For  the  calculations  reported  here 
we  have  chosen  an  Isothermal  equation  of  state  for  the  electrons 
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pp  I 

y ■ 1,  and  assumed  Chat  the  ions  are  cold  T.  » 
e i 

The  fluid  equations  are  in  Eulerian  form  and  are  solved  explicitly. 
The  densities  and  velocities  are  advanced  in  time  using  a modified 
Lax-Wendroft  two-step  procedure.2  In  the  spirit  of  Roberts  and  Weiss,3 
we  use  a fourth  order  algorithm  for  the  derivatives  of  the  fluxes  and 
a small  Courant  parameter,  At/Ax  >0.1.  The  algorithm  is  second  order 
in  time  and  fourth  order  in  space.  Poisson's  equation  is  solved 
utilizing  the  FFT.4 

Pereira5  has  also  solved  these  equations  numerically  in  his 
analysis  of  Langmuir  solitons.  His  objective  was  also  to  compliment 
work  involving  the  nonlinear  Schrodinger  equation.  He  observed  the 
propagation  of  these  solitons,  their  collisions,  and  their  etiology 
in  the  oscillating- two- stream  instability. 

We  initialize  the  ions  as  a travelling  ion  acoustic  sinusoidal 
potential  wave  with  a flow  velocity  equal  to  the  ion  acoustic  velocity 
in  the  direction  opposite  to  that  of  the  wave.  Disregarding  non- 
linearities,  the  result  is  a time  invariant  ion  density  modulation, 

n^  ■ 1 ■+»  6n^  sin(ksx).  (6) 

I The  ion  wave  has  a wavelength  of  bO 0 Debye  lengths  (i.e.  kg/kjj  = 2n/k00) 

and  therefore  there  will  be  a low  frequency  electron  density  fluctuation 
approximately  equal  to  this.  In  addition,  the  electrons  are  given  a 
shorter  wavelength  high  frequency  fluctuation.  For  example.  Fig.  1(a) 
shows  the  initial  electric  field  which  results  when  the  ion  wave  has 


an  amplitude  of  twenty  percent,  6n^  ■ 0.2,  and  Fourier  mode  19  is 

excited  as  a right  traveling  electrostatic  electron  plasma  wave  with 

a 0.1$  amplitude,  6n  » 0.001.  (The  electric  field  is  defined  as 

W/5x  and  is  therefore  scaled  to  V^n  T . ) For  reference  the 

o e 

discrete  Fourier  transform  of  this  wave  is  also  shown.  The  only 
nonzero  components  are  the  first,  the  phonon,  and  the  nineteenth,  the 
plasmon. 

Figure  1(b)  shows  configuration  space  and  k-space  at  time 
tu>p  - 108.  The  plasmon  has  evolved  to  have  shorter  wavelengths  in 
the  center  of  configuration  space.  In  k-spece  the  pUonon  is  unchanged, 
but  the  plasmon  is  now  represented  by  a broad  spectrum  of  modes.  If 
the  wave  is  launched  to  travel  in  the  opposite  direction,  the  mode 
spectrum  at  this  time  is  the  same,  but  the  localization  is  centered 
around  the  origin  in  confirmation  space. 

We  were  interested  in  whether  or  not  we  were  observing  a diffusion 
process  in  k- space  so  we  recorded  the  spread  of  the  Fourier  wave- 
numbers  of  the  plasmon.  This  is  defined 

Ak  - V<  k*  - < k >->  (7) 

where 

< kn  > » E (kQ|Ek|2)/£  |eJ2  , (8) 

lc  1c 


with  the  sum  extending  over  the  plasmon  wavenumbers. 


The  time  history  o£  this  spread  for  the  case  we  have  been  con- 
sidering is  shown  In  Fig.  2.  For  a scattering  process  which  can  be 


represented  by  a diffusion  equation  in  k-space,  the  square  of  the 
spread  (Ak)2,  the  variance,  will  be  proportional  to  time.  In  contrast 
Fig.  2 shows  Ak  to  be  approximately  proportional  to  time.  Furthermore, 
the  diffusion  theory  of  Davidson  and  Goldman6  indicates  that  the 
p lasmon-phonon  scattering  should  proceed  as  (6n^)2.  The  process  we 
are  observing,  however,  evolves  linearly  in  Sn^.  For  example.  Fig.  1(b) 
is  also  the  result  for  6n^  >0.1  (half  the  phonon  amplitude)  at  time 
tujp  =*  21 6 (twice  the  time)  with  the  other  parameters  unchanged. 

The  results  of  this  computation  show  that  the  high  frequency 

Fourier  modes  have  approximately  the  same  frequency  and  thus  maintain 

a constant  phase  relation  during  the  localization  process.  This  is 

0 

neglecting  times  when  the  amplitudes  change  sign,  then  there  are 
phase  jumps  of  magnitude  equal  to  tt.  With  this  in  mind  it  is  a 
straightforward  matter  to  obtain  these  results  analytically. 

The  electron  fluid  equations,  Eqs.  (1)  and  (2),  may  be  linearized 
and  combined  to  yield  the  following  equation, 


a2n 

— + (1  + 6^  sin  kgx) 


( ***•) 

\e  - Ye  -EF-/ 


(9) 


We  assume  that  the  initial  plasmon  distribution  is  sharply  peaked  in 

k-space  so  that  the  last  term  may  be  written  bzne/bx2  - - 

The  resulting  equation  may  be  solved  directly.  Using  Poisson's 
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equation  the  high  frequency  electric  field  is  given  by 
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6n 

E(x,t)  » E (x)  exp  [i(l  + sin  k x)uj  t]  (10) 

o da  po 

where  id  is  the  frequency  of  the  initial  electric  field  E (x), 
po  o ' 

u)2  » a)2  ( i + v k 2/kjr).  (We  have  assumeu  that  6n.  is  small  so 

po  p e e D i 

that  VI  + 6n^  *•  1 + 6n^/2.)  If  the  initial  condition  is  a right 
traveling  plasmon,  the  solution  i3 

6n 

E(x,t)  - Eq  sin  [kgx  - (1  + — ^ — sin  kgx)u>pot]  (11) 

We  can  define  a specially  dependent  effective  wavenumber  ke^ 
as  the  partial  of  the  argument  of  the  sine  with  respect  to  displacement. 
This  gives 


eff 
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k - k _ 

e s 2 
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-s-  (cos  k x)  u)  t. 

s po 


(12) 


This  wavenumber  increases  (the  plasmon  localizes)  near  the  center 

of  our  configuration  space  where  k x * tt.  This  is  the  situation  of 

s 

Fig.  1(b).  Similar  analysis  indicates  that  a left  traveling  plasmon 
will  localize  near  the  origin  as  we  observe. 

Returning  to  Eqs.  (7),  (8)  and  (10),  using  the  Bessel  generating 
function,  and  an  identity  in  Watson,7  it  is  straightforward  to  show 
that 


7 


2 


'13) 


(Ak)2 


[&k(e 


o)]2 


*(■ 


6n.k 
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For  Che  initially  monochromatic  case  this  result  may  be  obtained 
directly  from  Eq.  (12)  by  averaging  over  space.  In  Fig.  2 this  line 
has  been  added  as  the  dashed  line.  The  spread  In  the  Fourier  wave- 
numbers  initially  increases  according  to  this  equation. 

A batter  understanding  of  the  assumption  made  in  this  analysis 
is  obtained  by  substituting  our  solution  for  the  high  frequency 
density  fluctuation  into  Eq.  (9).  We  find  terms  proportional  to 
powers  of  time  which  do  not  cancel  since  the  electron  wave  does  not 
remain  monochromatic.  The  lowest  order  term  in  the  discrepancy  has 
a magnitude  equal  to  Veke($n^  k 

effect  of  the  density  modulation  and  neglected  this  term  which  is 
proportional  to  the  spread  in  k-space.  We  have  made  the  following 
assumption 


u) 

s po 


t)/kj. 


We  have  retained  the 


Veke(Ak)/kD  < 6ni/2  * UM 

The  left  hand-side  of  this  inequality  may  be  written  as 

Ak  6u)(k  * kg)/6k  which  is  a frequency  spread  due  to  a wavenumber 

spread.  When  this  spread  is  small  relative  to  that  caused  by  the 

density  modulation  then  our  analytic  solution  may  be  applied. 

We  wish  to  point  out  that  Eq.  (12)  predicts  another  interesting 

phenomenon.  At  the  time  to)  * 2k  /(k  6n.  ),  the  effective  wavenumber 

po  e s i 
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ke^c  wtll  go  to  zero  at  the  origin.  Henceforth,  the  plasmon  will  be 
trapped;  i.e.,  it  will  be  reflected  so  as  to  be  confined  to  certain 
regions  of  configuration  space.  In  our  simulations  we  have  observed 
this  but  our  understanding  here  is  limited.  Therefore,  another  limit 
to  this  theory  is  that  Ak  < kg. 

Our  essential  assumption  is  written  as  the  inequality  of  Eq.  (1^). 
This  led  to  a solution,  Eq.  (10),  which  in  turn  indicated  that  the 
spread  of  the  Fourier  wavenumbers  is  monotonically  increasing  in 
time,  Eq.  (15).  Our  theory  breaks  down  after  a period  of  time,  and 
is  therefore  valid  only  in  indicating  the  initial  response  of  the 
system. 

Figures  1(c)  and  lfd)  show  the  electric  field  and  its  Fourier 
spectrum  for  a smaller  ion  wave,  6n^  = 0.05,  somewhat  later  in  time. 

In  Fig.  1(c)  the  high  frequency  field  was  again  initialized  as  a 
right  traveling  plasmon  with  an  amplitude  5ng  = 0.001.  The  wave 
collapse  was  again  initially  near  the  center  of  our  configuration 
space,  but  at  ta^  33  250  it  began  slowly  moving  to  the  right.  This 
is  evidence  of  a balancing  of  the  last  two  terms  of  Eq.  (9).  The 
electron  fluid  of  Fig.  1(d)  was  initialized  with  modes  16-22  excited 
as  standing  waves,  all  with  an  amplitude  of  6ne  = 0.001.  Again  we 
observe  a balancing  of  the  final  terms  of  Eq.  (9).  In  this  case, 
with  several  Fourier  components  traveling  in  either  direction,  there 
is  clearly  an  envelope  structure  to  the  waves  in  configuration  space. 
Its  shape  and  position  correlate  with  that  of  the  ions  (|e|2  °=  - Sn^). 
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In  this  note  we  have  reported  the  results  of  our  numerical 
analysis  of  the  one-dimension  two-fluid  plasma  equations.  Our 
computations  are  formulated  as  an  initial-value  problem  and  the  initial 
behavior  of  our  system  of  Langmuir  oscillations  seems  to  contrast  with 
well-accepted  notions.  In  particular,  we  observe  the  Langmuir  waves 
to  collapse  or  localize  into  regions  of  strong  density  gradients  and 
not  into  density  depressions.  This  is  understood  by  noting  that  our 
particular  choice  of  initial  conditions  has  a narrow  spectrum  in  k- 
space.  As  such, variation  in  the  electron  pressure  term  may  be 
neglected  in  solving  for  the  initial  behavior.  If  this  is  done  the 
fluid  equations  have  a simple  analytic  solution  which  agrees  with  the 
results  of  our  simulations.  The  solution  indicates  that  the  width  of 
the  plasmon  spectrum  in  k-space  will  increase  in  time.  Later  in  time, 
after  our  analytic  solution  becomes  invalid,  the  position  of  the 
localization  moves  into  the  density  depression.  The  system  evolves 
into  a state  in  which  the  density  gradients  are  balanced  by  the 
electrostatic  force  and  the  force  of  spatially  varying  high  frequency 
electric  field.  These  later  states  exhibit  behavior  which  agree  with 
the  well-established  theoretical  notions. 
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(b)  twp  = 108 

Fig.  1 — The  electric  field  and  its  Fourier  transform,  (a)  An  initial  condition: 
nj  = 0.2  sin  ksx,  k,  ■ 2jrkD/400,  ne  = 0.001  sin  (kex— cupt),  ke  = 19kg.  (b)  This 
is  later  in  the  same  calculation,  toop  = 108.  The  collapse  of  the  Langmuir  rad- 
iation is  evident. 

12 


(d)twp*760 

Fig.  1 (Continued)  — The  electric  field  and  its  Fourier  transform,  (c)  In  this 
calculation  the  ion  wave  amplitude  is  reduced  to  n}  = 0.05;  the  electron  fluid 
has  the  same  initial  condition.  This  is  at  fu>p  * 710;  the  collapse  has  shifted 
to  the  right,  (d)  Again  the  ion  wave  has  an  amplitude  of  5%,  but  now  the 
electrons  were  initialized  with  a spectrum  of  standing  waves,  ne  * 0.001  2^ 
sin  (2k,x— 0g).  (0g  are  randomly  chosen  phase  angles.)  This  is  at  tu>p  * 

760;  the  field  clearly  exhibits  an  envelope  structure. 


